Implications of percolation of colour 
strings on multiplicities, correlations and 
the transverse momentum 

M.A.Braun 
High-energy physics department 
S.Petersburg University, 198904 S.Petersburg, Russia 

and C.Pajares 
Departamento de Ffsica de Partfculas, 
Universidade de Santiago de Compostela, 
15706-Santiago de Compostela, Spain 

July 1999 

Percolation, multiplicities, (p^) and correlations. 



Abstract. 

In the colour string model the impact of string percolation on multiplicities, (pj) 
and their long range (forward-backward) correlations is studied. It is assumed that dif- 
ferent string overlaps produce the observed hadrons independently. The multiplicities 
are shown to be damped by a simple factor which follows from the percolation theory. 
The (p 2 T ) rise at the same rate as multiplicities fall. A clear signature of the percolation 
phase transition is found to be a behaviour of the forward-backward correlations for 
intensive quantities, such as (p%) or its inverse, which can be detected in the relativistic 
heavy ion collider. 
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1 Introduction 

Multiparticle production at high energies is currently (and successfully) described in 
terms of colour strings stretched between the projectile and target [1-6]. Hadronization 
of these strings produces the observed hadrons. Colour strings may be viewed upon 
as (small) areas in the transverse space filled with colour field created by the colliding 
partons. Creation of particles goes via emission of qq pairs in this field. With growing 
energy and/or atomic number of colliding particles, the number of strings grows. Once 
strings have a certain nonzero dimension in the transverse space they start to overlap 
forming clusters, very much like disks in the 2-dimensional percolation theory. The 
geometrical behaviour of strings in the transverse plane then follows that of percolating 
discs. In particular at a certain critical string density a macroscopic cluster appears 
(infinite in the thermodynamic limit), which marks the percolation phase transition [7 
-9]. 

The percolation theory governs the geometrical pattern of the string clustering. 
Its observable implications however require introduction of some dynamics to describe 
string interaction, that is, the behavior of a cluster formed by several overlapping 
strings. 

One can study different possibilities. 

A most naive attitude is to assume that nothing happens as strings overlap, in other 
words, they continue to emit particles independently without being affected by their 
overlapping neighbours. This is a scenario of non-interacting strings, which closely 
corresponds to original calculations in the colour strings approach, oriented at com- 
paratively small energies (and numbers of strings). This scenario however contradicts 
the idea that strings are areas of the transversal space filled with colour field and thus 
with energy, since in the overlapping areas the energy should have grown. 

In another limiting case one may assume that a cluster of several overlapping strings 
behaves as a single string with an appropriately higher colour field (a string of higher 
colour, or a "colour rope" [10]). This fusion scenario was proposed by the authors and 
later realized as a Monte-Carlo algorithm nearly a decade ago [11,12]. It predicts low- 
ering of total multiplicities and forward-backward correlations (FBC) and also strange 
baryon enhancement, in a reasonable agreement with the known experimental trends. 

However both discussed scenarios are obviously of a limiting sort. In a typical 
situation strings only partially overlap and there is no reason to expect them to fuse 
into a single stringy object, especially if the overlap is small. The transversal space 
occupied by a cluster of overlapping strings splits into a number of areas in which 
different number of string overlap, including areas where no overlapping takes place. 
In each such area colour fields coming from the overlapping strings will add together. 
As a result the total cluster area is split in domains with different colour field strength. 
As a first approximation, neglecting the interaction at the domain frontiers, one may 
assume that emission of qq pairs in the domains proceeds independently, governed by 
the field strength ("the string tension") in a given domain. This picture implies that 
clustering of strings actually leads to their proliferation, rather than fusion, since each 
particular overlap may be considered as a separate string. Evidently newly formed 
strings differ not only in their colours but also in their transverse areas. 

As a simple example consider a cluster of two partially overlapping strings (Fig. 1). 
One distinguishes three different regions: regions 1 and 2 where no overlapping takes 
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place and the colour field remains the same as in a single string, and the overlap region 
3 with colour fields of both strings summed. In our picture particle production will 
proceed independently from these three areas, that is, from three different "strings" 
corresponding to areas 1,2 and 3. In this sense string interaction has split two strings 
into three of different colour, area and form in the transverse space. 

We stress that these dynamical assumptions are rather independent of the geomet- 
rical picture of clusterization. In particular, in each of the scenarios discussed above, 
at a certain string density there occurs the percolation phase transition. However its 
experimental signatures crucially depend on the dynamical contents of string interac- 
tion. With no interaction, clustering does not change physical observables, so that 
the geometric percolation will not be felt at all. With the interaction between strings 
turned on, clustering (and percolation) lead to well observable implications. 

In this note we shall study these implications for simplest observables, such as 
multiplicities, average transverse momenta and FBC in the realistic scenario discussed 
above, which corresponds to independent particle production from different overlap 
domains. Our choice of observables is dictated by the possibility to relate them directly 
to the domain properties, without introducing any more assumptions. 

2 Multiplicity and (p^) for overlapping strings 

As stated in the Introduction, the central dynamical problem is to find how the observ- 
ables change when several strings form a cluster partially overlapping . In the overlap 
areas the colour fields of individual strings are summed together. It is more convenient 
to sum the charges which generate the colour field of overlapping strings. 

Let only two strings, of areas oo in the transverse space each, partially overlap in 
the area S 2 (region 3 in Fig. 1), so that Si = <jq — S 2 is the area in each string not 
overlapping with the other. In the following it will be called the overlap area of one 
string. A natural assumption seems to be that the average colour density £ of the 
string in the transverse plane is a constant 

i = Qo/vo, (i) 

where Qo is a color of the string. For partially overlapping strings the colour in each 
of the two non-overlapping areas will then be 

Q l =ZS 1 = Q (S l /a ). (2) 

The colour in the overlap area Q 2 will be a vector sum of the two overlapping colours 
C,S 2 . In this summation the total colour squared should be conserved [10]. Indeed, 
Q\ = (Qov + Qou) 2 where Q ov and are the two vector colours in the overlap area. 
Since the colours in the two strings may generally be oriented in an arbitrary manner 
respective to one another, the average of Q i;Qo„ is zero. Then Q\ = Ql v + Q' ov 2 , which 
leads to 

Q 2 = V2£S 2 = V2Q (S 2 /a ). (3) 

One observes that, due to its vector nature, the colour in the overlap is less than the 
sum of the two overlapping colours. This phenomenon was first mentioned in [10] for 
the so-called colour ropes. 
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Two simplest observables, the multiplicity // and the average transverse momentum 
squared (p?) are directly related to the field strength in the string and thus to its 
generating colour. In fact they are both proportional to the colour [10,13]. Thus, 
assuming independent emission from the three regions 1,2 and 3 in Fig. 1 we get for 
the multiplicity 

fi/fi = 2(S 1 /a ) + V2(S 2 /a ), (4) 

where /io is a multiplicity for a single string. To find (p^) one has to divide the total 
transverse momentum squared of all observed particles by the total multiplicity. In 
this way for our cluster of two strings we obtain 

/_2w/_2v _ 2(S 1 /a ) + V2V2(S 2 /a ) _ 2 

\Pt//\Pt/o 2 (S 1 /a Q ) + V2(S 2 /a Q ) 2{S 1 /a Q ) + ^ (S 2 /a ) ' KJ 

where (pj)o is the average transverse momentum squared for a single string and we 
have used the evident property 2 Si + 2S2 = 2o"o in the second equality. 

Generalizing to any number N of overlapping strings we find the total multiplicity 

as 

/V^Ev^^M)), (6) 

i 

where the sum goes over all individual overlaps i of rii strings having areas Simi- 
larly for the (p\) we obtain 

In the second equality we again used an evident identity J2i n iS^ — Na Q . Note that 
(6) and (7) imply a simple relation between the multiplicity and transverse momentum 

= N, (8) 

Mo (Pt)o 

which evidently has a meaning of conservation of the total transverse momentum pro- 
duced. 

Eqs. (6) and (7) do not look easy to apply. To calculate the sums over i one seems 
to have to identify all individual overlaps of any number of strings with their areas. 
For a large number of strings the latter may have very complicated forms and their 
analysis presents great calculational difficulties. However one immediately recognizes 
that such individual tracking of overlaps is not at all necessary. One can combine all 
terms with a given number of overlapping strings rii = n into a single term, which sums 
all such overlaps into a total area of exactly n overlapping strings S n . Then one finds 
instead of (6) and (7) 

N 

/VVo = J2Vn(S n /a ) (9) 

71=1 

and 

( P 2 t)/(p 2 t)o = ^ N Lo, y (10) 

E n =l V^(S n /(To) 
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In contrast to individual overlap areas the total ones S n can be easily calculated 
(see Appendix). Let the projections of the strings onto the transverse space be dis- 
tributed uniformly in the interaction area S with a density p. Introduce a dimensionless 
parameter 

rj = pa = Na /S. (11) 

The "thermodynamic limit" corresponds to taking the number of the strings N — > oo 
keeping r] fixed. In this limit one readily finds that the distribution of the overlaps of 
n strings is Poissonian with a mean value r?: 

Pn ~ S ~ n\ 6 ■ [iZ) 
From (9) we then find that the multiplicity is damped due to overlapping by a factor 

< 13 » 

where the average is taken over the Poissonian distribution (12). 

The behaviour of F(rj) is shown in Fig. 2. It smoothly goes down from unity at 
rj = to values around 0.5 at rj = 4 falling as for larger 77's. According to (10) 

the inverse of F shows the rise of the (pt)- Note that a crude estimate of F{rj) can be 
done from the overall compression of the string area due to overlapping. The fraction 
of the total area occupied by the strings according to (12) (see also [14]) is given by 

]r Pn = l-e-". (14) 

n=l 

The compression is given by (14) divided by 77. According to our picture the multiplicity 
is damped by the square root of the compression factor, so that the damping factor is 



1 - e- r i 

nn) = ] ]—^-. (is) 

For all the seeming crudeness of this estimate, (15) is very close to the exact result, as 
shown in Fig. 2 by a dashed curve. 



3 Percolation 

Percolation is a purely classical mechanism. Overlapping strings form clusters. At some 
critical value of the parameter rj a phase transition of the 2nd order occurs: a cluster 
appears which extends over the whole surface (an infinite cluster in the thermodynamic 
limit). The critical value of 77 is found to be i] c ~ 1.12 — 1.20 [15]. Below the phase 
transition point, for 77 < ?7 C , there is no infinite cluster, Above the transition point, at 
i] > i] c an infinite cluster appears with a probability 

Poo = 0(v-Vc)(v-Vcf- (16) 

The critical exponent j3 can be calculated from Monte-Carlo simulations. However the 
universality of critical behaviour, that is, its independence of the percolating substrate, 
allows to borrow its value from lattice percolation, where (3 = 5/36. 
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Cluster configuration can be characterized by the occupation numbers (v n ), that is 
average numbers of clusters made of n strings. Their behaviour at all values of rj and 
n is not known. From scaling considerations in the vicinity of the phase transition it 
has been found [24] 

(u n ) =n- T F(n a (ri-r] c )), \q - r] c \ « 1, n » 1, (17) 

where r = 187/91 and a = 36/91 and the function F(z) is finite at z = and falls 
off exponentially for \z\ — > oo. Eq. (17) is of limited value, since near rj — r\ c the 
bulk of the contribution is still supplied by low values of n, for which (14) is not valid. 
However from (17) one can find non-analytic parts of other quantities of interest at the 
transition point. In particular, one finds a singular part of the total number of clusters 
M = Yi, v n as A(M) = c\i] — i] c \ 8 / 3 . This singularity is quite weak: not only (M) itself 
but also its two first derivatives in rj stay continuous at rj = r\ c and only the third blows 
up as \r\ — r] c \~ l l z . So one should not expect that the percolation phase transition will 
be clearly reflected in some peculiar behaviour of standard observables. 

Indeed we observed that neither the total multiplicity nor (p\) show any irregularity 
in the vicinity of the phase transition, that is, at rj around unity. This is not surprizing 
since both quantities reflect the overlap structure rather than the cluster one. The 
connectedness property implied in the latter has no effect on these global observables. 

It is remarkable, however, that the fluctuations of these observables carry some 
information about the phase transition. The dispersion of the multiplicity due to 
overlapping and clustering can easily be calculated in the thermodynamic limit (see 
Appendix). The result is shown in Fig. 3. The dispersion shows a clear maximum 
around rj = 1 (in fact at rj ~ 0.7). So some information of the percolation phenomenon 
is passed to the total multiplicity, in spite of the fact that it basically does not feel the 
connectedness properties of the formed clusters. Of course, due to relation (10), the 
dispersion of (pj) has a similar behaviour. 

We have to warn against a simplistic interpretation of this result. The dispersion 
shown in Fig. 3 is only part of the total one, which besides includes contributions from 
the fluctuations inside the strings and also in their number. Below we shall discuss the 
relevance and magnitude of these extra contributions. 

An intriguing question is a relation between the percolation and formation of the 
quark-gluon plasma. Formally these phenomena are different. Percolation is related to 
the connectedness property of the strings. The (cold) quark-gluon plasma formation is 
related to the density of the produced particles (or, equivalently, the density of their 
transverse energy ). However in practice percolation and plasma formation go together. 
In fact, the transverse energy density inside a single string seems to be sufficient for the 
plasma formation. Percolation makes the total area occupied by strings comparable 
to the total interaction area, thus, creating, a sizable area with energy densities above 
the plasma formation threshold. 

Let us make some crude estimates. Comparison with the observed multiplicity 
densities in pp(p) collisons at present energies fix the number of produced (charged) 
particles per string per unit rapidity at approximately unity. Taking the average energy 
of each particle as 0.4 GeV (which is certainly a lower bound), formation length in the 
Bjorken formula [16] as 1 fm and the string transverse radius as 0.2 fm [7] we get the 
3-dimensional transverse energy density inside the string as ~ 3GeV/fm 3 . The plasma 
threshold is currently estimated to be at lGeV/im 3 . So it is tempting to say that the 
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plasma already exists inside strings. This however has little physical sense because of 
a very small area occupied by a string. One can speak of a plasma only when the total 
area occupied by a cluster of strings reaches a sizable fraction of the total interaction 
area. In Fig. 4 we show this fraction for a maximal cluster as a function of r\ calculated 
by Monte -Carlo simulations in a system of 50 strings. It grows wit rj and the fastest 
growth occurs precisely in the region of the percolation phase transition: as rj grows 
from 0.8 to 1.2 the fraction grows from 0.3 to 0.6. With a string cluster occupying 
more than half of the interacting area, one can safely speak of a plasma formed in that 
area. 



4 Dispersions and forward- backward correlations 

As shown in the previous section a definite signal for the percolation phase transition 
comes from the multiplicity (or ipj)) fluctuations due to string clustering. However 
these fluctuations are not the only ones. An important contribution also comes from the 
fluctuations of the multiplicity inside the strings, or rather inside the individual overlaps 
of strings, which in our picture play the role of independent particle emitters. To find 
the total dispersion it is convenient to use a formalism of the generating functions. Let 
the total probability to observe n produced particles be V(n). In our picture it is given 
by a convolution of the probability for a given overlap configuration P(C) with the 
probabilties for particle production from all individual overlaps: 

V(n) = J2P(C) ]T p 1 (n 1 )...p M (n M )6 ntEni . (18) 

C m,...n M 

Here C is a "configuration", characterized by the total number of the overlaps M 
and their individual properties: area and number of overlapping strings. We pass to 
generating functions 

*(z) = 5> B P(n), h(z) = 5> B Pi(n), (19) 

n n 

for which (18) transforms into 

Hz) = Y,P(cmMz)- (20) 

c i 

The overall averages (with the probability V) are given by 

<n> = &(z) z=1 , (n(n-l)) = $"(z) z=1 . (21) 
The averages inside the i-th overlap are likewise given by 

m = <P'{z) z=1 , mim - l) = <P"{z) z=1 . (22) 

Using this formulas one readily finds the total dispersion of multiplicity as a sum of 
two terms 

D 2 = D 2 C + Dl. (23) 
Here D 2 C is a dispersion due to fluctuation in configurations 

D 2 c = H P{C) E^™-(E p (C)E^) 2 - (24) 

C ik c i 
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It is calculated according to the distribution V assuming that the numbers of particles 
produced in individual overlaps are fixed to be their averages. It is this part of the 
dispesion, which is calculated in the Appendix and presented in Fig.3. The part D 2 n 
is a part of the dispersion due to fluctuations inside the individual overlaps: 

^ = EmEK-(^)' (25) 

c i 

Its calculation is extremely difficult even in by Monte Carlo simulations, since it requires 
identification of all individual overlaps and knowledge of their areas. To simplify we 
assume that the distribution of produced particles from any overlap is Poissonian, so 
that nf — {jh) 2 = n[. Then we evidently find 

Dl = ^ (26) 

From the experimental data we estimate /xq — 11 for unit rapidity interval. Com- 
parison of Fig. 2 and Fig. 3 then shows that for unit rapidity interval the internal 
dispersion D 2 n is roughly 40 times greater than the "percolation dispersion" D 2 C at 
i] = 0.7 corresponding to the maximum for the latter. At other 77 the ratio D 2 n /Dfj is 
still greater. With the growth of the rapidity interval this ratio diminishes proportion- 
ally. However it is obvious that one should not expect to clearly see the percolation 
effects directly in the observed multiplicities. 

A better signal for the percolation comes from the FBC, which, as we shall see, 
distinguish between the percolation and intrinsic dispersions (in fact they depend on 
their ratio). The FBC are described by the dependence of the average multiplicity in 
the backward hemisphere (//b) on the event multiplicity in the forward hemisphere fi F - 
The data can be fitted by a linear expression [1] 

(/i B ) = a + bfj, F , (27) 
where a and b are given by expectation values 



(/4> - M 2 



(28) 



6 - /-,2\ /„ \2 • ( 29 ) 



In absence of FBC (hf^b) = (i i f)(i ji 'b) and one obtains a = (fxp) and 6 = 0. So the 
strength of the correlations is given by the coefficient b. 

To calculate the necessary averages we introduce the probability V(F, B) to produce 
F{B) particles in the forward (backward) hemispheres. Similar to (18) it is given by a 
convolution 

V(F,B) = £P(C) £ Hp^BJS^Sz^ (30) 

C* Fi,Bi i 

which transforms into a relation between the generating functions 

$(z F z B ) = Y J P{C)\{<P{z F ,z B ). (31) 

C i 
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Instead of (21) we now have 



8zf I i V dz% I , \ dzpdzs I i 

(32) 

formulas similar to the first two for B, and similar formulas for the averages over the 
distributions pi inside the overlap i with the generating function <fii(zF, zb)- Using this 
formulas we find that the dispersions again split into the percolacion (C) and internal 
(in) parts. 

D 2 F = D 2 FC + D 2 Fin , D FB = (FB) - (F) (B) = D FB C + D 2 FBin . (33) 
They are given by 

d%c = E no E Fi ^ - (E PIP) (34) 

C i.k C i 

D 2 fb,c = E P (C) E Fi * ~ E P (C) E Fi E p (c) E ^ (35) 

C* i,fc C* i C i 

D F , m = J2P(C)E(F?-(W), (36) 

C* i 

D 2 FB, m = E p (^) E(^ - ^ (37) 

C i 

The usual assumption is that there are no FB correlations for particle production 
from a single emitter (overap % in our case). Then D 2 FBin = 0. Also from n = F + B 
one can relate these dispersions with the overall ones D 2 C and D 2 n : 

Dl,c = Dfb,c = (V4)^, D 2 F . n = (l/2)D 2 n . (38) 
From (29) one then finds 

1 D 2 

So, the FBC parameter b indeed measures the ratio of internal to percolation dispersions 
squared. If we assume that the distribution of particles produced from an individual 
overlap is Poissonian then (39) transforms into 



1 n , o l± 



Both n and D 2 C grow linearly with the number of strings N, so in the thermodynamic 
limit b does not depend on N. However it depends on the multiplicity of a single string 
Ho, since \x is proportional to fiQ and D 2 C is proportional to /j/q. As mentioned, from 
the experimental data we find that /io — 1.1 y where y is the rapidity interval of 
the produced particles. So the right-hand side of (40) falls with the growth of the 
rapidity window. It also has a minimum at i] ~ 0.7 corresponding to the maximum 
D 2 C . From this we see that the FBC parameter b also has a maximum near rj ~ 0.7 
whose magnitude grows with the rapidity interval. In Fig. 5 we show the behaviour of 
b as a function of i] for two rapidity intervals y = 5 and y = 9 
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One can also study the FBC for (p?). Due to relation (8), for fixed N they are 
uniquely determined by the fluctuations in //. It is convenient to choose the inverse 
(Pt) as an observable: 

t = ^p¥T = n- (41) 

Then it is obvious that the parameter b for r is the same as for the multiplicity, since 
both percolation and internal dispersions for r are obtained from those for //, divided by 
N. However this equivalence holds only for fixed number of strings N. In the realistic 
situation N fluctuates. The multiplicity is an extensive observable and carries these 
fluctuations directly. They result large and smear out nearly all traces of percolation. 
In contrast, r is an intensive observable and so feels the fluctuation in N only indirectly, 
through the fluctuations in the parameter 77. As a result, the percolation effects for the 
FBC in r are clearly visible, as we shall see in the next section. 



5 Realistic hadronic and nuclear collisions 

For realistic hadronic and nuclear collisions the total number of strings iV grows with 
energy and fluctuates. For pp collisions the initial number of colour strings N created 
at a given cm. energy yfs can be taken from the well-known calculations of [1,2]. 
For pA and AB collisions this number should by multiplied by the effective number of 
collisions C. At asymptotic energies for minimum bias hA and AB collisions the AGK 
rules predict 

C A = Aa/a A , C AB = ABa/a AB , (42) 

where a (a A) a AB ) is the inelastic pp (pA, AB) cross-section. For hA collisions at fixed 
impact parameter b 

C A = AaT A (b)/(l-exp(-AaT AB )), (43) 

where T A is the standard nuclear profile function (normalized to unity). The same 
relation holds for AB collisions in the optical approximation with 

A - AB, T A - T AB (b) = J d 2 cT A {c)T B {b - c). (44) 

The transverse interaction area for pp and hA collsions obviously is of the order a. 
For AB collisions it depends on the geometry of the collisions. We restrict ourselves 
to central (b = 0) collisions of identical nuclei, when evidently the interaction area is 
a AA . The parameter 77 is then calculated according to (11). We have taken the string 
area a as 

a = na 2 , a = 0.2 fm. (45) 

in accordance with arguments presented in [7] . To relate the multiplicity to the exper- 
imental data we normalized fi to the observable central (charged) multiplicity per unit 
rapidity in pp collisions at low energies. This fixes fiQ for unit rapidity to be 1.1 The 
found values of fx = (dn ch /dy) y=0 are presented in the Table for central S-S and Pb-PB 
collisions (4d column) together with the corresponding values of c.m.energy y/s, r\ and 
number of strings iV (1st to 3d columns, respectively). To compare, we present the 
values of \x found without fusion and percolation (independent strings picture) in the 
7th column. One should have in mind that the asymptotic formulas for the number 
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of collisions (42) - (44) used to determine the number of strings are not valid at com- 
paratively low energies due to restrictions imposed by energy conservation. As stated 
in [1], at lower energies the number of collisions can be roughly obtained multiplying 
(42) - (44) by a factor 1/2. Correspondingly, our values for the multiplicity at two 
lower energies have to be corrected for energy conservation by a factor of this order. 
We preferred not to make this correction, since, in any case, it cannot be determined 
with any degree of rigour, and the effects we are considering are essential only at high 
enough energies, where (42)- (44) are hopefully valid. 

As we observe, percolation considerably damps the multiplicity at high energies: at 
t/s = 7000 GeV for central Pb-Pb collision from nearly 8000 down to approximately 
3000. This effect was predicted in our earlier papers on the string fusion [11,12] and is 
now reproduced in various models [17]. 

As to the FBC, we studied them both for the multiplicity and the inverse (jpj) 
(in fact for the observable r, defined in (41)). We have taken one half of the total 
rapidity available for the relevant rapidity window. The "external" dispersion in the 
denominator of (39) has to include also the fluctuations in the number of strings. We 
have assumed the overall distribution in N to be Poissonian so that AN/N = l/\/N. 
As mentioned, for the multiplicity, due to its extensive character, the dispersion in N 
simply adds to Dc- It is large and absolutely dominates all other contributions. So it 
is no surprise that the coefficient b for the multiplicity with percolation (5th columns 
in the Table) only slightly differs from the one without percolation (8th column). 

In contrast, for r the fluctuations in N only enter via the dependence of rj on 
N. They result smalller than the percolation contribution at the maximum of the 
latter but fall with rj more slowly, so that they dominate at large rj. As a result, the 
characteristic peaked form of the parameter b found at fixed N (Fig. 5) disppears 
and b results steadily growing with 77 (6th column in the Tables). However one should 
note that in the independent string picture the parameter b for r is zero, since the 
dependence on N is then completely absent. So experimental study of the FBC for 
inverse (pf,) seems to be a promising way to observe signatures of string fusion and 
percolation. 

6 Conclusions 

In this study we have analyzed the impact of fusion and percolation of colour strings 
on global observables, such as multiplicities and (pf,) . To do this, a certain dynamical 
assumption has been made. The strings have been assumed to decay into the observed 
hadrons independently in each overlap. 

On the qualitative level the results are best seen in the idealized case of a fixed 
number of strings N >> 1 (equivalent to plasma studies in the thermodynamic limit). 
A clear consequence of fusion and percolation is damping of the multiplicities, well 
described by a damping factor (15), which follows from the percolation theory. An 
unexpected but potentially important result is that the parameter b of the FBC shows 
a clear maximum at the percolation point. Such a maximum would be natural in 
fluctuations, say, of the cluster sizes, where it is to be expected as a signature of the 
percolation phase transition. However multiplicities do not seem to feel directly the 
cluster structure, so that the appearance of the maximum in their fluctuations is a new 
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result. 

In the realistic case of nuclear collisions, where N fluctuates, the predicted multi- 
plicities repeat the pattern observed at fixed N and are damped by the same factor 
(15). At the LHC energies it reduces the multiplicities by more than two times. How- 
ever the parameter b for multiplicity results completely dominated by fluctuations in 
the number of strings and so only slightly different from the one in the independent 
string picture. A better (intensive) observable to see the impact of percolation is (p\) 
or its inverse, for which the contribution of the fluctuations in N is drastically reduced. 
For (pt)~ 1 we predict sizable values of b to be contrasted with b = in the independent 
string picture. Observation of nonzero values of b for (p 2 -) -1 would therefore be a clear 
signature of string fusion and percolation. 
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8 Appendix. Multiplicities and their dispersion 

Let us recall the geometry of our percolation picture. Discs of radius a and area 
o"o = 7ra 2 are homogeneously distributed in the total area S. It is assumed that centers 
of the discs are inside the unit circle of area So = it so that S = 7r(l + a) 2 . The 
disc density is p = N/S and the percolation parameter is 77 = pa = Na /S. In the 
thermodynamic limit — > 00, so that at fixed rj the radius of the discs goes to zero. 
For fixed r\ 

• =yf-i)\ («•> 



so that at large N a ~ 1/yN and a ~ 1/N. Since the discs are distributed homoge- 
neously, the probability that their centers are at points rj, % = 1, ...N inside the unit 
circle is independent of r\ and is given by 

Pin) = S» N . (47) 

Let us take a configuration which corresponds to the disc centers at points r\. Then 
the overlap area of exactly n discs is given by the integral 

S n (n,...r N )= fd 2 r ]T f[6(a-\r-r lk \) {[ 0(\v - v ik \ - a). (48) 

{h,..in}C{ii,...i N } k=l k=n+l 

The average of Sn will be given by a multiple integral over r\ with the probability (47): 

N 



where 



1 r JV r 
<^> = ^v/ U<P ri S n ( ri ,...r N ) = C% d 2 rF n (r)(l-F(r)f- n , (49) 

F(r) = (l/S )( d 2 ri e(a-\r (50) 

J Sn 
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The function SoF(r) gives an area occupied by a circle C of radius a with a center 
at r which is inside the unit circle So. If r < 1 — a then C is always inside So so that 

F(r) = a /S , < r < 1 - a. (51) 

However for r > 1 — a a part of C turns out to be outside the unit circle, and 

F(r) = cr(r)/S , 1 - a < r < 1 + a, (52) 

where a(r) < <jq is the overlap of the two discs C and So- 

Generally an overlap of two circles of radii r\ and r 2 with a distance r between 
their centers is given by 

a(r 1} r 2 ,r) = (l/2)r^(a 1 - sin ax) + (1/2)^(^2 - sina^), (53) 

where 



cos(ai/2) = 
cos(o;2/2) = 




The function a(r) in (52) is just a(l,a,r). 

Formulas (49) -(54) allow to calculate numerically the average (S n ) for any finite 
value of iV without much difficulty 

In the thermodynamic limit, iV — > oo with r\ being fixed, the calculation of (S n ) 
becomes trivial. Indeed then one can neglect the part of integration in r with r > 1 — a 
altogether, with an error ~ l/a ~ l/y/N. With the same precision one then finds 



(S n ) = SC n N (a G /S) n (l - a /S) N ~ n , (55) 

where we have put So — S The physically relevant values of n remain finite as iV — > oo. 
So we can approximately take 

C n N = N n /n\, (1 - a /S) N ~ n = exp(-Na / S). (56) 



We then find that in the thermodynamic limit the distribution of overlaps in n is 
Poissonian with the mean value given by i] (Eq. (13)). 

Calculation of the multiplicity dispersion requires knowledge of the average of its 
square. With the centers of the discs at r±, ...rjv it has the form 

/i 2 (n, ...r N ) = (l/a 2 )CT V^S n (r u ...r N ))\ (57) 

n 

where Sn(ti, ...rjy) is given by (48). Taking the average over the discs centers positions 
we now come to a double integral in r and r' 

1 r 1 r N 

= Z2 E d 2 rd 2 r'^ Y[d 2 n 

°0 m,n Jb °0 Jb o i=l 

n N 

E II%-|r-rd) n 0(\v-v ik \-a) 

{h,..i„}C{l,...N} k=l k=n+l 
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m N 

E 11%-lr-rJ) II 0{\v-v n \-a). (58) 

{jl,-Jm}C{l,.. .«}«=! Z=m+1 

This complicated expression, however, continues to be factorized in all rj. Let us 
assume for the moment that n > m. We can always rename the variables to have 
in the first sum r 1? ...r n as variables in the first product of ^-functions (and as a result, 
the rest N — n variables r ra+1 , ..rjy go into the second product). We shall then have 
Cpj terms with the identical first sum. Now let p variables in the first product of 6- 
functions in it coincide with p variables from the set r±, ...r n . We have such terms, 
which will evidenly all have the same dependence on the mentioned p variables. The 
left m — p variables from the first product of ^-functions in the second sum do not 
coincide with any variables r±, ...r n , so they are chosen fron variables r n+ i, ...tn- We 
shall have C^Zn various terms of this sort. Thus the overall symmetry factor turns out 
to be a multibynomial coefficient 

r"i nV s^rn-p _ _ (~<p,n-p,m-p /r Q \ 

p\yn — py.ym — p)\[n — n — m + p)\ 

It multiplies the result of integration over all r iy % = 1,...N, which has the form (at 
fixed r and r') 

P (r, r>) X n - p (r, r') X m - p (r', r )C N ~ n ~ m+p (r, r'), (60) 

where 

<t>( r y) = (l/S ) / 9(a-\r-r 1 \)9(a-\r'-r 1 \), 



I So 

X (r, r') = (1/S ) I 6(a - |r - nD^lr' - n| - a) = F(r) - 0(r, /), 

J So 

((r,r) = (1/S ) [ e(\r-T 1 \-a)e(\T f -T 1 \-a) = l-F(r)-F(r')+<f>(r,r'), (61) 

J So 

with F(r) defined before by (50). Of course one should sum over all possible values of 
p. 

Combining all the terms, we find the expression for the average square of multiplicity 

as 

= K E ^(n + p)(n + p) 



°0 n. 



m,p 



CT* I d 2 rd 2 r'r(r, r') X n (r, r') X m (r', r)^""" m " p (r, r'). (62) 
J s 

This expression is exact and may serve as a basis for the calculation of the average 
square of the multiplicity at finite N. However the new function becomes very 
complicated when both variables r and r' are greater than 1 — a (it is then given by 
the overlap area of three circles and we do not know any simple analytic expression for 
it). 

For this reason rather than analyze the general expression (62) for finite iV we shall 
immediately take the thermodynamic limit iV — > oo. We are in fact interested in the 
dispersion, not in the average square of multiplicity. It is important, since the leading 
terms in N cancel in the dispersion. So we shall study the difference 



D 2 = ^) - (/i) 2 



(63) 
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in the limit N — > oo, r] finite. As we shall see, although both terms in the right-hand 
side of of (63) behave as iV 2 separately, their difference grows only as N. 

Separating from (62) the term with p = and combinig it with the second term 
on the right-hand side of (63) we present the total dispersion squared as a sum of two 
terms 

D 2 = Dj + D\, (64) 



where 



D\ = - 2 Y,^I d 2 rd 2 r'[C^ m x n (r,r') X m (r',r)C N - n - m (r,r')- 

On „ ™ J S 



a n,m 

C%C%F n (r)F M (r')(l - F(r)) N - n (l - F(r')) N - m ] (65) 

and D\ is given by (62) with a restriction p > 1. 

Function S (f>(r, r') gives the overlapping area of three circles: the unit circle So 
and two circles C and C of radii a with centers at r and r' . It is evidently zero if 
R = |r — r'| > 2a since for such R circles C and C do not overlap. Because of this, in 
the part D 2 , which has at least one factor 0, the integration in r and r' is restricted to 
the domain R < 2a, whereas in the part D\ the integration covers the whole range of 
values of R. 

We begin with the part D\. Splitting the integration region in r and r' in two parts, 
R > 2a and R < 2a we make use of (61) and present the first part as 

D 2 X1 = ^V v& / d 2 rd 2 r'F n (r)F m (r') 
°o n,m J R> 2a 

[C% m (l - F(r) - F(r')f~ n - m - C%C%(1 - F(r)f~ n (l - F(r')f~ m }. (66) 

One notices immediately that the in the thermodynamic limit the leading terms in the 
integrand (independent of N) cancel and only terms of the order 1/N remain. The 
integration over r and r' provides a factor oc iV 2 so that the total contribution results 
oc N. To find this term we use that at large N, up to terms of the order 1/iV 2 , 

_ N n ( n(n-l) \ 

CiV ~ ^ 2AT- J' 

rm ,n _ N n+m ( {m + n){m + n-l) \ 
Cn -^r\ L 2N )> 

(1 _ F )N-n = e ~N F{l +nF _ NF 2/2), 

where we have taken into account that F has order 1/N. Then we obtain 

D 2 = -1 V ^-N n+m [ d 2 rd 2 r'F n (r)F m (r') exp(-N(F(r) + F(r')) 
w nlml jR>2a 

(mF(r) + nF{r') - NF(r)F(r') - nm/N). (67) 

The integrand is now explicitly of the order 1/N, so that we can change the integration 
region to r,r' < 1 — a, since the difference in the area will be of the order a 2 ~ 1/N, 
which results in the overall difference of the order 1/N 2 and can safely be neglected. 
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In the region r,r' < 1 — a both F(r) and F(r) are constants, given by (51). Integration 
over r and r' gives an overall factor S 2 , so that in the end we get 

D 11 /N — —e '} — — -y/nmynln + m) — ri — nm\ = . (68) 

In the last expression the averages are to be taken over the Poissonian distribution 
with the mean value 77. The part D\ x is thus negative (and results comparatively small 
for all values of rj. 

In the second part of D\ the integration goes over a small region R < 2a, of the 
order a 2 ~ 1/N, so that one can retain only the leading terms in the integrand. Passing 
to the integration over r and R one observes that at fixed r < 1 — a the integration 
region over R covers the whole region R < 2a. For r > 1 — a the integration region in 
R becomes much more complicated, determined by the condition r' — |r + R| < 1 + a. 
However the contribution from the latter region will be of the order a 3 ~ 1/N\/N, 
since apart from a factor oc a 2 from the integration over R, a factor oc a appears due to 
integration over r > 1 — a. So, up to terms of the relative order l/y/N, we can neglect 
the contribution from the region r > 1 — a. 

If r < 1 — a the circle C is completely inside the unit circle Sq. Then its intersection 
with the circle C is also automatically inside S . Therefore function S <f>(r,r) in this 
region is simply given by the overlap of the circles C and C, that is a(a, a, R) defined 
by (53). We thus find 

<t>(r,r') = ^\(R), (69) 

where 

\(R) = (l/n)(a - sin a), a = 2 arccos(i?/2). (70) 

Putting this into the expression for D\ 2 and retaining the leading terms in the limit 
iV — > 00 we obtain 

D\ 2 /N = -e~ 2r > £ y^. v n+m f 2 RdR[(l - \{R)) n+m e^ - 1]. (71) 

V n,m n - m - JO 

This expression can be easily evaluated numerically. It is relatively large and also 
negative. 

We finally come to the part D\. The integration in r,r' goes over R < 2a, so that 
we can apply the same approximations as made in calculating D\ 2 . We find 



D 2 2 /N = -e- 2r > Yl E ^ n+ f^ + -r] n+m+ P f RdRX>(R)(l - \(R)) n+m e^ R l 
V n,m p =i n.m.p. Jo 

(72) 

This part is evidently positive. Its numerical evaluation shows that it nearly cancels 
the large negative contributions from D\ (in fact four digits are cancelled typically). 
So the numerical calculation of the dispersion requires some care. 
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10 Figure captions 

Fig. 1. Projections of two overlapping strings onto the transverse plane. 
Fig. 2. Damping of the multiplicity as a function of rj. 

Fig. 3. Percolation dispersion squared of the multiplicity per string in units /1q as 
a function of r\. 

Fig. 4. Fraction of the total interaction area covered by the maximal cluster as a 
function of 77. 

Fig. 5. Parameter b of the FBC for different rapidity intervals as a function of 77 
for fixed N. 



11 Table captions 

Columns show from left to right: cm. energy per nucleon *Js\ parameter 77 (Eq. (1)); 
number of initially produced strings N; central charged multiplicity per unit rapidity 
fi and its FBC parameter B(/j,), the FBC parameter for the inverse (pj<), b(r), and the 
multiplicities and 6(/ii) in the independent string model (subscript 0). 



Table 

S-S scattering (b = 0) 
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N 






b(r) 




bo 


19.4 


0.33 


99 


101 


0.67 


0.04 


111 


0.65 


62.5 


0.47 


143 


140 


0.74 


0.08 


160 


0.72 


200 


0.65 


198 


185 
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0.13 


221 
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546 
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255 


228 
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284 


0.80 


1800 


1.10 


336 


283 


0.85 


0.24 


374 


0.82 


7000 


1.47 


448 


349 


0.87 


0.30 


499 


0.85 



Pb-Pb scattering (b = 0) 





V 


N 


A* 


%) 


b(r) 




bo 


19.4 


1.45 


1530 


1200 


0.69 


0.13 


1700 


0.65 


62.5 


2.09 


2200 


1530 


0.75 


0.19 


2450 


0.72 


200.0 


2.88 


3040 


1870 


0.78 


0.25 


3390 


0.77 


546.0 


3.72 


3920 


2170 


0.80 


0.28 


4370 


0.80 


1800.0 


4.90 


5170 


2530 


0.81 


0.31 


5760 


0.82 


7000.0 


6.53 


6890 


2940 


0.82 


0.33 


7680 


0.85 



13 Figures 




Figure 1: Projections of two overlapping strings onto the transverse plane. 




Figure 2: Damping of the multiplicity as a function of 77. 
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Figure 3: Percolation dispersion squared of the multiplicity per string in units 
function of 77. 
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Figure 5: Parameter b of the FBC for different rapidity intervals as a function of r\ for 
fixed N. 



